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摘 要 : 准确 评价 地 下 水 对 河水 的 补给 量 是 流域 水 资源 量 管理 和 合理 利用 的 基础 。 在 马鞍 河流 域 
下 游 采 集 不 同 季节 地 表 水 和 地 下 水 样品 75 组 ,利用 Cl (EC) FDO 同位 素 多 方法 联合 
评价 ,识别 了 地 下 水 补给 河水 的 位 置 、 补 给 量 及 其 季节 变化 。 结 果 表 明 : 马 鞍 河 水 EC 和 Cl 质量 浓 
度 沿 着 流向 均 呈 降低 趋势 ,3D 和 S O 值 沿 流向 减 小 。 雨 季 EC 和 Cl 质量 浓度 最 低 ,8D m S 0 M 
最 高 。 地 下 水 各 组 分 浓度 均 低 于 河水 ,时 空 变 化 不 明显 。 地 下 水 单 宽 排泄 量 存在 时 空 变 异 , 上 段 和 


段 为 地 下 水 强 排泄 区 ,中 段 地 下 水 排泄 较 弱 ,不 同 季节 地 下 水 排泄 量 占 总 排泄 量 的 72.20% ~ 


下 
95.07% 。 雨 季 地 下 水 单 宽 排 泄 量 显著 降低 ,河水 中 基 流 比例 由 雨季 前 期 的 68.89% 降 至 29.43% 。 
整体 上 ,地 下 水 补给 河水 季节 变化 明显 ,而 空间 变化 规律 较为 稳定 。 研 究 成 果 有 利于 深入 认识 河水 


和 地 下 水 的 相互 作用 机 制 ,并 为 当地 水 资源 利用 提供 基础 依据 。 
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地 下 水 补给 河水 是 地 表 水 和 地 下 水 相互 作用 关 
系 的 一 种 中 ,在 干旱 半 干 旱 区 普遍 存在 。 由 于 地 下 
水 水 质 水 温 相 对 稳定 ,其 补给 对 河流 水 文系 统 和 生 
态 系统 具有 重要 调节 意义 1。 准确 获取 地 下 水 补 
给 河水 的 位 置 和 补给 量 , 是 水 资源 调查 评价 的 基础 ， 
也 是 流域 水 资源 合理 配置 的 前 提 ”= 。 由 于 水 体 的 
化 学 组 分 和 同位 素 特征 记录 了 水 体形 成 .循环 和 演 
化 的 历史 中 成 为 定量 评价 地 下 水 补给 河水 的 有 效 
手段 , 张 兵 等 5 -5 学 者 以 水 化 学 和 氧 氧 稳 定 同 位 素 
为 示 踪 指标 ,分 别 研究 了 松花 江 、 巴 音 河 海流 图 河 、 
昭 苏 一 特 克 斯 例 地 流域 水 循环 特征 和 地 表 水 与 地 下 
水 转化 关系 。 

然而 ,地 表 水 和 地 下 水 相互 作用 受到 多 种 复杂 
因素 影响 ,单一 方法 评价 结果 存在 较 大 的 不 确定 
性 ” ,需要 开展 多 方法 联合 示 踪 并 相互 验证 以 提高 
评价 精度 "1 。 另 外 ,已 有 研究 多 在 单一 时 段 评 价 ， 
对 地 下 水 补给 河水 的 动态 变化 关注 较 少 。 实 际 上 ， 
这 种 补 排 关系 受 水 动力 场 控制 ,降雨 导致 的 河水 和 
地 下 水 流 场 季节 性 变化 必然 会 影响 两 者 交换 量 …) ， 
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分 析 地 下 水 补给 河水 的 季节 变化 有 利于 深入 认识 河 
水 和 地 下 水 作用 关系 。 

马鞍 河 位 于 黄土 高 原 腹地 ,是 甘肃 省 庆 阳 市 境 
内 的 主要 地 表 水 系 , 在 当地 水 资源 格局 中 具有 重要 
地 位 。 苏 小 四 、 侯 广 才 等 ”学 者 在 该 地 区 的 研究 
表明 浅 层 地 下 水 补给 马 莲 河 ,但 均 以 定性 分 析 为 主 。 
本 次 工作 在 前 人 研究 成 果 的 基础 上 ,综合 利用 水 体 
Cl -电导 率 (EC) 和 O 等 指标 ,基于 质量 平衡 原理 
定量 评价 了 地 下 水 补给 河水 的 时 空 变化 , 旨 在 为 水 
资源 合理 利用 提供 依据 ,并 为 同类 研究 提供 参考 。 


1 研究 区 概况 


马 莲 河 发 源 于 陕西 省 白 于 山 一 带 , 干 流 流 经 甘 
肃 省 庆 阳 市 ,于 宁县 汇 信 泾 河 。 总 长 度 375 km , 流 
域 面积 1.92 x 10° km ,多 年 平均 径流 量 4.51 x 10° 
m ,河流 宽度 5.3 ~48.1 m, A EPER] FY ZH 
加 。 研 究 区 为 庆 城 县 以 南 的 流域 下 游 ,地 处 上 暧 温带 
半 湿 润 气候 区 ,年 平均 降水 量 570 mm ,水 面 蒸发 量 
1565 mm。 地 表 乡 为 黄土 覆盖 ,地 形 切割 强烈 ,地 表 
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图 1 研究 区 取样 点 分 布 和 水 文 地 质 剖 面 图 


Fig. 1 Sampling locations and hydrogeological profile of the study area 


水 系 发 育 ,流量 较 大 的 支流 包括 合 水 川 和 九龙 河 。 
研究 区 黄土 地 下 水 零星 分 布 , 仅 在 面积 较 大 的 
黄土 堪 区 具 供 水 意义 。 第 四 系 松散 岩 类 孔隙 水 呈 狭 
长 状 分 布 在 河谷 阶地 ,含水 层 厚 度 多 小 于 5 m, 水 量 
极为 贫乏 。 白 琴 系 地 层 沉积 厚度 巨大 , 陆 相 雄 悄 岩 
TRARRE A , 赋 存 规模 地 下 水 ,成 为 本 区 分 布 最 广 
泛 的 含水 岩 组 (图 1) 。 总 体 上 ,可 将 其 分 为 上 部 环 
河 组 含水 层 ( Kih) 和 下 部 洛 河 组 含水 层 (Ki1) 。 环 
河 组 含水 介质 为 泥 质 砂岩 和 细 砂 岩 , 厚 度 约 300 m, 
地 下 水 水 位 埋 深 从 黄土 震 区 的 200 m 至 河谷 区 的 
5 m。 受 新 构造 运动 的 影响 , 马 莲 河 持续 下 切 ,切割 
下 伏 的 环 河 组 含水 层 。 地 下 水 从 两 侧 分 水 岭 流向 河 
谷 , 通 过 河床 与 河水 产生 水 力 联系 。 洛 河 组 含水 层 
含水 介质 为 粗 砂 岩 , 属 深 循 环 地 下 水 ,交替 缓慢 。 


2 研究 方法 


2.1 取样 和 测试 

2016 年 分 别 于 雨季 前 (4 月 中 旬 ) .雨季 (8 月 中 
旬 ) 和 雨季 后 (11 月 中 旬 ), 沿 着 马 莲 河 下 游 庆 城 至 
入 河口 约 100 km 的 河 段 采集 河水 .地 下 水 和 支 沟 水 
样品 ,用 于 测试 EC .CI 和 和 氧 氧 稳定 同位 素 。 河 水 
取样 点 选择 在 易于 取样 的 水 文 站 和 公路 桥 等 位 置 ， 
采样 间距 5 ~10 km, 河水 取 自 河流 中 部 河床 0.3 m 
以 上 部 位 。 

沿岸 采集 地 下 水 样品 ,采样 位 置 尽量 靠近 河流 
以 反映 补给 特征 , 受 城镇 分 布 格局 影响 ,取样 的 地 下 


水 井 均 位 于 河流 左岸 。 支 沟 水 样品 采 自 板 桥 水 文 站 
和 新 宁 水 文 站 ,分 别 代 表 下 游 主要 支 沟 合 水 川 和 九 
龙 河 。 每 一 期 次 采集 河水 样品 15 组 .地 下 水 样品 8 
组 和 支 沟 水 样品 2 组 ,共计 3 个 期 次 75 组 。 每 个 时 
段 采样 时 均 避 开 降雨 事件 ,并 在 1 ~2 d 内 完成 。 

EC 参数 利用 便携 式 多 参数 水 质 分 析 仪 (HACH 
HQ40D ) 现 场 测试 ,Cl “0 和 D 测试 由 中 国 地 质 
调查 局 水 文 地 质 环境 地 质 调查 中 心 实验 室 完成 ,其 
中 Cl 采用 ICS2000 型 离子 色谱 仪 测定 (误差 
<2.0% ) ,8D 和 5"0 采用 波长 扫描 光 强 误 荡 光谱 
技术 测试 ,精度 分 别 为 0. 10%o 和 0. 03%o, 分 析 结 
采用 VSMOW 标准 。 马 莲 河 干流 上 庆 阳 、 宁 县 、 雨 落 
坪 等 水 文 站 流量 数据 来 自 黄 河水 利 委 员 会 西峰 水 文 
水 资源 勘测 局 , 支 沟 板 桥 .新 宁 水 文 站 流量 数据 来 自 
庆 阳 水 文 局 。 

2.2 模型 

取样 时 段 马 莲 河 沿途 工农 业 取 水 较 少 ,予以 忽 
略 不 计 。 对 任 一 取样 段 ,河水 流量 均衡 要 素 包 括 上 
游 来 水 ,沿途 地 下 水 汇流 、 降 雨 产 流 、 茸 发 和 支 沟 补 
给 。 其 中 ,降雨 产 流 较 为 复杂 ,由 于 本 次 研究 中 未 开 
展 相关 工作 , 产 流 的 计算 参考 已 有 研究 成 果 。 卢 龙 
彬 等 '” 在 黄土 高 原 北部 研究 发 现 , 径 流 系 数 主要 受 
平均 降雨 强度 影响 ,并 存在 如 下 统计 学 关系 : 

a = 0.047R* + 0.175R (1) 


式 中 :a 为 径流 系数 ;R 为 平均 降雨 强度 [mm ， 
(5 min) ~'), 


tó wy 


| 


— 


h=0, #¥%0.1 利用 公式 (4) 利用 公式 (5) 
Sy EL. I+ 计算 Cr Cn HEU Um 


图 2 7 值 计算 流程 图 


Fig.2 Flow chart for calculating the value of groundwater inflow rates 


3 结果 与 分 析 


得 出 取样 段 降雨 产 流 为 : 
Qp = aPS (2) 


式 中 :0; 为 降雨 产 流 (m d7); P 为 日 降雨 量 
(mm + d7') ;5 为 取样 段 内 扣除 支 沟 流域 面积 后 的 
产 流 面积 (mm ) 。 
由 此 , 任 一 取样 段 水 量 平衡 方程 如 下 : 
Q,= Q +I, +Q,+ Qp- EW, (3) 
式 中 :0 OO, 为 取样 区 间 起 点 ,终点 河流 流量 、 
沟 补 给 量 (m + do!) ;7 为 地 下 水 单 宽 排泄 量 ( 研究 
区 地 下 水 补给 河水 ,1>0,m . do!) ;x 为 取样 区 间 
距离 (m) ;五 为 河流 蒸发 强度 (m+ d”) ;下 为 河流 
DBT TE it FE Cm.) 。 
对 河流 中 某 一 化 学 或 同位 素 组 分 ,假设 在 河水 
中 混合 均匀 , 若 忽 略 沿途 水 岩 作 用 及 同位 素 分 馏 等 
因素 影响 ,质量 平衡 方程 如 下 : 
Q, C= OC, + L Cg, +QpCp +Q, C, (4) 


IP : Ci, Ca, Cun Cp. C, 为 取样 起 点 河水 、 终 点 河 
水 `. 地 下 水 .降雨 和 支 沟 水 中 某 组 分 质量 浓度 。 

根据 公式 (4) ,对 于 每 一 取样 区 间 ,C, 和 了 存在 
定量 关系 。 利 用 Cl” EC 和 "0 3 种 指标 ,建立 取样 
点 河水 各 组 分 模拟 值 和 实测 值 的 误差 函数 ,定义 
如 下 : 


(5) 


UD -3(“"] 


ô, 


3.1 EC 和 Cl 的 时 空 变化 

为 分 析 各 组 分 质量 浓度 的 空间 差异 ,将 马 莲 河 
下 游 分 为 3 段 , 庆 城 县 至 合 水 县 为 上 段 , 合 水 县 至 宁 
县 为 中 段 ,宁县 以 下 为 下 段 。 年 内 降雨 分 布 和 水 文 
情势 具有 明显 的 季节 变化 , 据 宁 县 气象 站 和 雨 落 坪 
水 文 站 监测 数据 ,1 ~5 月 流域 内 降雨 较 少 ,马鞍 河 
流量 相对 稳定 。6 ~9 月 降雨 增多 ,暴雨 频 发 , 马 莲 
河流 量 明 显 增高 。10 ~ 12 月 ,降雨 量 减少 ,河水 流 
量 逐 渐 降 低 并 趋 于 稳定 。 

根据 测试 结果 ( 表 1) ,2016 年 马 莲 河水 EC 值 
1705.40 ~ 3 123. 60 ps .cm ,Cl 质量 浓度 
387. 20 ~707.42 mg - Lo! ,远大 于 长 江 、 黄 河 及 珠江 
流域 "~" 。 因 为 马 莲 河上 游 河水 受 高 矿 化 度 地 下 
水 补给 ,具有 较 高 的 背景 值 。 受 此 影响 ,下 游 河 
水 EC 值 和 Cl 质量 浓度 也 较 高 。 各 组 分 质量 浓度 
季节 变化 明显 (图 2) ,基本 规律 是 雨季 后 期 最 高 .两 
季 前 期 次 之 、 雨 季 最 低 。4 月 雨季 前 期 河水 EC 值 
1 753.20 ~2 869.36 ps + cm ,均值 2 329.41 ps - 
em”! ,Cl~ 质量 浓度 401.30 ~ 639.39 mg + 工 -均值 
485.23 mg- L, 4 8 A mÆ, EC 值 1 705. 40 ~ 
2 642.00 us + cm ,均值 2 244.47 us .cm ,Cl 
质量 浓度 387. 20 ~ 531.38 mg ' 工 一 ,均值 469. 29 
mg ' 工 - 。 可 以 看 出 ,相对 雨季 前 期 ,雨季 河水 EC 
值 和 Cl 质量 浓度 均 有 所 降低 ,这 可 以 归结 为 雨季 


式 中 :VU(7) 为 误差 函数 ;C;、C 为 取样 点 某 组 分 质量 
浓度 模拟 值 和 实测 值 (i 取 1、2.3 ,分 别 代表 3 种 指 
ER) ;6; 为 某 组 分 质量 浓度 标准 差 ,用 于 消除 不 同 指 
标量 纲 影 响 。 

对 每 一 取样 段 , 以 河水 C SEC 和 ”0 值 为 拟 合 
目标 ,以 地 下 水 单 宽 排泄 量 7 为 拟 合 参数 ,利用 公式 
(5) 计 算 误 差 函 数 。 设 置 7 初 始 值 为 0 , 步 长 为 0.1， 
运行 模型 直至 误差 函数 最 小 (图 2) ,获取 雨季 前 期 、 


河水 流量 增高 导致 的 稀释 作用 。 到 11 月 雨季 后 期 ， 
河水 EC 值 和 Cl 质量 浓度 有 所 增加 ,均值 分 别 为 
2 444.87 us + cm "和 536.81 mg + L , 高 出 雨季 前 
期 。 可 能 是 因 雨 季 后 期 随 着 降雨 的 减少 ,河水 离子 
组 分 受到 蒸发 浓缩 影响 而 增高 。 支 沟 水 EC 值 
1 031.60 ~ 2 984. 53 us .cm ,Cl 质量 浓度 
412.30 ~599.26 mg . L ,同样 是 8 月 雨季 最 低 。 
地 下 水 EC 值 409. 58 ~1 696.84 us + cm ,Cl 质 


雨季 和 雨季 后 期 3 个 时 段 地 下 水 单 宽 排泄 量 的 空间 


量 浓度 115. 92 ~- 255. 96 mg L71, 均 低 于 河水 ,不 


分 布 。 


同时 段 各 取样 点 EC 值 和 Cl- 质量 浓度 呈 波 动 状 , 季 
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表 1 不 同时 段 河水 和 地 下 水 EC、Cl 和 和 氢 氧 同位 素 统计 


Tab.1 Average Cl” , EC and hydrogen and oxygen isotopes of river water and groundwater during three periods 


河水 地 下 水 
时 段 组 分 = = = = 
最 小 值 最 大 值 均值 最 小 值 最 大 值 均值 

4 月 EC /hs:cm-! 1 753.20 2 869. 36 2 329.41 410.90 1 696. 84 836.33 

Cl-/mg + L7! 401.30 639. 39 485.23 130.52 243. 23 173.71 

580 /%o -9.96 -8.75 -9.37 -10.93 -9.26 -10.29 

8D /%o -74.55 -66.97 -71.55 -80.80 -71.51 -76.04 
8 月  EC/ps* cm"! 1 705.40 2 642.00 2 244.47 409.58 1 654. 30 820.21 

Cl-/mg+ L7! 387. 20 531.38 469. 29 122.31 216. 10 172.64 

80 /%o -7.45 -6.30 -6.83 -10.54 -9.11 -9.78 

8D /%o -61.96 -54.88 -58.40 -78.91 -68.28 -72.66 
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河水 和 地 下 水 Cl 质量 浓度 和 电导 率 EC 的 时 空 变化 


Fig.3 Spatial and seasonal variations of Cl” and EC in river water and groundwater 


节 变 化 不 明显 。 

从 图 3 可 以 看 出 , 沿 着 河水 流向 ,各 取样 时 段 河 
水 EC 值 和 Cl 质量 浓度 除了 在 中 段 略 有 波动 外 , 基 
本 上 呈 降 低 趋势 ,地 下 水 各 组 分 无 明显 空间 变化 规 
律 。 一 般 来 说 ,在 以 贿 石 风化 或 蒸发 结晶 作用 为 主 
的 地 区 , 沿 着 流向 河水 EC 值 和 Cl 质量 浓度 不 断 增 
加 ,而 马 莲 河 河水 却 显示 出 沿 着 流向 降低 的 相反 
趋势 。 考 虑 到 相对 于 河水 流量 ,沿途 支流 补给 较 少 ， 
河水 EC 值 和 C1 质量 浓度 不 断 降低 的 原因 应 是 沿 
途 持续 接受 低 EC 和 Cl 质量 浓度 的 地 下 水 汇流 
补给 。 


3.2 和 氢 氧 同位 素 的 时 空 变化 

据 测试 结果 ( 表 1) ,4 月 河水 5D 和 5"0 值 分 
别 为 -74.55%o ~ -66.97%o 、-9.96%o ~ — 8. 75%o, 
8 月 分 别 为 - 61. 96%o ~ - 54. 88%o 、- 7. 45%o ~ 
-6.30%o,11 月 分 别 为 - 78. 25%o ~ - 68. 85%o、 
-10.41%o ~ -9.05%o。 可 以 看 出 ,河水 氧 氧 稳定 同 
位 素 特征 具有 显著 的 季节 变化 ,53D 和 a O 值 在 雨 
季 明 显 增高 ,雨季 后 期 降低 ,这 和 河水 的 补给 来 源 及 
补给 源 水 的 氢 氧 同位 素 特 征 有 关 。 雨 季 河 水 接受 大 
量 大 气 降雨 补给 ,河水 同位 素 特征 主要 受降 雨 同位 
素 特征 有 影响。 研究 区 地 处 黄土 高 原 南 部 ,雨季 水 汽 
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主要 来 源 于 东部 太平 洋 和 六 盘山 局 地 水 循环 ,已 
有 研究 发 现 黄土 高 原 地 区 大 气 降水 同位 素 在 夏季 富 
集 "1。 受 此 影响 ,雨季 河水 也 具有 较 高 的 SD 和 
80 值 。 地 下 水 4 月 8D 和 880 值 分 别 为 -80. 8%o 
~ -70. 51%0, — 10. 93%o ~ -9.26%o,8 月 分 别 为 
~78.91%0 ~ -68.28%o 、- 10. 54%o ~ -9.11%o,11 
月 分 别 为 - 79. 80%o ~ -76. 47%o 、- 11. 19%o ~ 
-10. 54%o。 地 下 水 8D 和 ò" O 值 低 于 河水 ,由 于 补 
给 来 源 较 为 稳定 , 氧 氧 同位 素 特征 季节 变化 不 明显 。 
图 4 显示 , 沿 着 河水 流向 ,不 同时 段 河 水 SD 和 5"0 
值 均 呈 现 出 下 降 趋 势 ,同样 为 沿途 接受 了 低 BD 和 
3 o 值 的 地 下 水 补给 所 致 。 沿 河水 流向 ,地 下 水 氧 
氧 同位 素 组 分 无 明显 变化 规律 。 

由 于 缺少 降雨 同位 素数 据 , 采 用 前 人 建立 的 当 
地 大 气 降 雨 线 ,方程 为 SD =7. 035380 -1.057 。 
从 氧 氧 同位 素 关系 图 (图 5) 可 以 看 出 ,尽管 不 同时 
期 地 下 水 均 沿 着 当地 大 气 降雨 线 分 布 ,但 多 位 于 降 
雨 线 左 下 端 , 且 各 时 段 地 下 水 8D 和 8° O 值 均 低 于 


河水 ,指示 地 下 水 和 河水 可 能 具有 不 同 的 补给 来 源 。 
研究 区 地 处 黄土 高 原 , 较 多 研究 表明 ,由 于 厚 
层 包 气 带 的 影响 ,黄土 高 原 地 下 水 主要 来 自 较 早 时 
期 降水 补给 ,现代 降雨 补给 微弱 ,地 下 水 同位 素 特征 


存在 偏 负 现象 。 河 水 和 支 沟 水 均 位 于 大 气 降雨 线 下 
方 , 且 不 同 程 度 偏离 ,指示 河水 经 历 蒸发 作用 。 男 
外 ,雨季 河水 SD 和 5"0 值 高 于 雨季 前 期 和 雨季 后 
期 ,这 是 因为 雨季 降雨 较为 集中 ,河水 对 大 气 降雨 响 
应 迅速 ,接受 了 5D 和 5"0 值 较 高 的 降雨 补给 ,造成 
同位 素 富 集 。 同 时 ,河流 为 连续 水 体 ,下游 同位 素 特 
征 受 到 上 游 影 响 。 马 莲 河 上 游 位 于 半 干 旱 气 候 区 ， 
蒸发 较为 强烈 ,导致 河水 重 同位 素 富 集 ,起 点 处 河水 
R01 样品 具有 最 高 的 SD 和 S O 值 即 为 优 证 。 
3.3 ”地 下 水 补给 河水 的 时 空 变化 

由 于 地 下 水 组 分 存在 空间 变异 ,上 段 .中 段 和 下 
段 地 下 水 端 元 分 别 取 各 段 采 集 的 地 下 水 样品 Cl 、 
EC 和 "0 均值 ,模型 所 需 各 参数 确定 如 下 : 马 莲 河 
初始 流量 Q, 和 支 沟 补给 量 分 别 取 自 庆 阳 \ 板 桥 、 新 
宁 水 文 站 取样 期 间 监 测 数据 .河流 蒸发 强度 取 4.90 
mm + d ,河流 断面 宽度 下 为 每 段 实 测 值 19. 62 ~ 
48.10 m, 河 水 Cl” .EC 和 "0 初始 值 分 别 为 R01 样 
品 测试 结果 。 

通过 对 比 R02 ~ R15 取样 点 各 指标 模拟 值 和 实 
测 值 判别 模型 评价 精度 ,图 6 显示 了 各 取样 点 Cl 、 
EC 和 "0 值 模拟 值 和 实测 值 拟 合 程度 。 可 以 看 出 ， 
各 样品 均 沿 着 1: 1 线 分 布 ,表明 拟 合 效果 较 好 。 其 
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图 4 河水 氢 氧 同位 素 的 时 空 变化 


Fig.4 Spatial and seasonal variations of hydrogen and oxygen isotopes in river water 
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图 5 不 同时 段 水 体 氧 氧 同 位 素 关 系 


Fig.5 Plots of 8D versus 8'*O during three periods 
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中 雨季 前 期 .雨季 和 雨季 后 期 Cl 质量 浓度 模拟 值 
和 实测 值 相对 误差 分 别 为 0.02% ~5.46% ,0.05% ~ 
7.35% 和 0.28% ~4.09% ,EC 模拟 值 和 实测 值 相 
对 误差 分 别 为 1.20% ~ 12.99% ,0.57% ~11.26% 
和 0.01% ~7. 11% ,0 模拟 值 和 实测 值 相对 误差 
分 别 为 0.06% ~2.63% ,0.06% ~3.78% 和 0.20% ~ 
3.80% ,各 时 段 大 多 数 样品 模拟 值 和 实测 值 相对 误 
差 低 于 5% 。 男 外 ,利用 马 莲 河 干流 宁县 、 雨 落 坪 水 
文 站 同期 监测 流量 和 模拟 流量 对 比 发 现 (图 7b ) , BY 
季 前 期 宁县 站 、 雨 落 坪 站 模拟 流量 相对 误差 分 别 为 
7.49% 和 1.85% ,雨季 相对 误差 分 别 为 3. 05% 和 
1.57% ,雨季 后 期 相对 误差 分 别 为 4. 91% 和 
14.59% 。 对 比 各 取样 点 Cl” LEC 和 "0 值 模拟 值 和 
实测 值 水 文 站 模拟 流量 和 监测 流量 ,显示 相对 误差 
较 小 ,表明 评价 结果 具有 可 信 性 。 这 也 说 明 多 方法 
联合 示 踪 可 作为 定量 评价 地 下 水 补给 河水 的 一 种 有 
效 手段 。 该 方法 基于 质量 守恒 原理 ,通过 沿 河 取样 
可 精准 识别 短 时 间 尺 度 地 下 水 排泄 位 置 和 单 宽 排泄 
量 , 在 水 文 监测 资料 较 少 的 流域 具有 应 用 优势 。 

根据 模拟 结果 (图 7) ,雨季 前 期 .雨季 和 雨季 后 
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王 雨 山 等 :基于 水 化 学 和 同位 素 评价 马 莲 河 下 游 地 下 水 补给 河水 的 时 空 变化 


8.4m +d! 0~12.9 m -d7', > Fie Bk B 
AREE 5lm + s7! 2.71m? + s7'\5. 10 
。 地 下 水 排泄 存在 空间 变异 ,上 段 (RO1 ~ 
Ra ee ee ee 
分 别 为 6.1 ~10.6m .di'!'1.5~5.3m © d', 
3.9 ~12.9 m? + d” ,排泄 量 分 别 占 地 下 水 总 排泄 量 
的 39% ,23.8% ,35.3% ; 中段 (R04 ~ R11) 段 单 宽 
排泄 量 较 低 ,该 段 河 流 长 度 占 下 游 河 流 总 长 度 的 
44.2% ,不 同时 段 地 下 水 排泄 量 仅 占 总 排泄 量 的 
4.9% ~13.8% ;下 段 (R11 ~ R14 ) 为 地 下 水 强 排 汇 
段 , 雨 季 前 期 .雨季 和 雨季 后 期 地 下 水 单 宽 排 汇 量 分 
别 为 7.1~13.0m +d! 1.3 ~8.4 m +d" 7.1 ~ 
11.1 m? + d”, 排 泄 量 分 别 占 地 下 水 总 排泄 量 的 
56.1% ,62.3% ,39.9% 。R14 ~ R15 段 地 下 水 排泄 
极为 微弱 ,不同 时段 单 宽 排 泄 量 0 ~3.9 m + do", 
各 时 段 模拟 河水 流量 及 基 流 比例 沿途 均 呈 增加 趋势 
(图 7b ,图 8) , 仅 在 地 下 水 弱 排泄 段 受 蒸发 影响 而 
略 有 降低 。 
地 下 水 排泄 同样 存在 季节 变化 ,雨季 地 下 水 单 
宽 排泄 量 显著 降低 (图 7a) ,雨季 前 期 和 后 期 变化 不 
大 ,图 8 显示 ,整个 下 游 河水 中 基 流 比例 在 雨季 前 
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图 6 各 指标 模拟 值 和 实测 值 对 比 


Fig.6 Comparison of modeling results and field measurements for Cl 


地 下 水 单 宽 排 汇 量 /m2* qd- 


沿 河 距 离 /km 
图 7 模拟 地 下 水 单 宽 排 泄 量 和 河水 流量 的 时 空 变化 


Fig.7 Spatial and seasonal variations of modeled groundwater inflow rates and river discharge 
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图 8 基 流 比例 的 时 空 变 化 
Fig.8 Spatial and seasonal variations of the ratios of 


baseflow in total river discharge 


期 和 后 期 分 别 为 68.9% 71.2% ,雨季 这 一 比例 将 
至 29.4% 。 尽 管 雨季 河水 流量 增加 明显 ,但 地 下 水 
总 排泄 量 降低 ,河水 流量 的 增加 主要 来 自 上 游 和 支 
沟 补 给 。 这 是 因为 相对 于 河流 对 降雨 的 迅速 响应 ， 
地 下 水 响应 存在 滞后 。 雨 季 河 水 水 位 增高 较 快 ,而 
地 下 水 水 位 则 相对 稳定 ,导致 含水 层 一 河流 之 间 水 
力 梯度 降低 ,在 其 他 条 件 不 变 情 况 下 ,地 下 水 对 河水 
的 补给 量 也 随 之 减少 。 人 研究 还 发 现 ,尽管 不 同时 段 
地 下 水 单 宽 排泄 量变 化 较 大 ,但 其 空间 变化 规律 较 
稳定 , 强 排泄 段位 置 未 发 生 改 变 。 这 是 由 于 地 下 水 
补给 河水 主要 由 含水 层 分 布 .地 下 水 动力 场 和 河床 
渗透 系数 等 因素 决定 ,而 这 些 因 素 随 时 间 变 化 较 小 。 

地 下 水 补给 河水 的 空间 变化 和 研究 区 地 质 水 
文 地 质 条 件 有 关 。 马 莲 河 下 游 上 段 和 下 段 均 穿 越 黄 
土 残 娘 ,含水 层 切割 强烈 ,地 下 水 水 力 梯 度 高 ,排泄 
量 大 。 中 段 马 莲 河 两 侧 地 形 较为 平坦 ,含水 层 切 割 
较 弱 ,地 下 水 水 力 梯度 低 ,排泄 量 小 。 另 一 方面 , 河 
床 岩 性 也 有 影响 ,在 砂岩 分 布 区 河床 渗透 系数 较 大 ， 
有 利于 地 下 水 补给 河水 ,而 在 泥 质 分 布 区 河床 渗透 
系数 较 低 ,地 下 水 和 河水 水 力 联系 较 弱 。 

地 表 水 和 地 下 水 转化 的 定量 评价 是 水 文学 领域 
的 研究 难点 ,理论 上 通过 断面 测 流 可 以 精确 评价 
河水 和 地 下 水 的 转化 过 程 ,但 在 实际 中 难以 实现 。 
本 次 研究 通过 利用 水 化 学 和 同位 素 多 方法 联合 示 
踪 , 评 价 了 马 莲 河流 域 地 下 水 补给 河水 的 时 空 变化 ， 
并 通过 和 实测 流量 .实测 河水 指标 对 比 ,证 明了 该 方 
法 具有 可 信和 性 。 但 限于 含水 层 一 河水 的 非 均 质 性 和 
模型 概 化 等 因素 影响 ,评价 结果 存在 不 确定 性 。 首 
先 , 地 下 水 端 元 的 代表 性 ,地 下 水 各 组 分 在 水 平 向 和 
重 向 上 存在 空间 变异 ,如 何 获取 排泄 区 的 地 下 水 样 


本 次 研究 中 以 各 段 采集 的 地 下 水 样品 C1” EC 


和 "0 均值 作为 上 段 ,中段 和 下 段 地 下 水 端 元 ,存在 
一 定 的 误差 。 其 次 ,河水 的 取样 密度 也 会 影响 评价 
结果 ,取样 密度 越 大 越 能 准确 刻画 地 下 水 补给 河水 
的 空间 变化 。 下 一 步 研究 中 ,应 根据 河水 水 化 学 和 
同位 素 指标 的 变化 特征 优化 取样 工作 。 


4 结论 

(1) BMk EC 值 1 705.4 ~3 123.6 ps - 
cm ,C1 质量 浓度 387.2 ~707.4 mg :上 L , 沿 着 流 
向 呈 降 低 趋 势 。 季 节 变 化 明显 ,雨季 后 期 最 高 .雨季 
前 期 其 次 雨季 最 低 。 地 下 水 EC 值 409. 6 ~ 
1 696.8 ps + cm !, Cl” 质量 浓度 115. 9 ~ 256. 0 
mg + L7! , 均 低 于 河水 ,季节 变化 不 明显 。 

(2) 不 同时 段 河 水 8D 和 8"0 值 沿 着 流向 呈现 
出 下 降 趋 势 ,雨季 河水 接受 SD 和 5"0 值 较 高 的 降 
雨 补给 , 重 同位 素 富 集 ,雨季 后 期 降低 。 地 下 水 8D 
和 6"0 值 低 于 河水 ,无 明显 季节 变化 规律 。 

(3) 地 下 水 单 宽 排 汇 量 存在 时 空 变化 ,上 段 和 
下 段 为 地 下 水 强 排泄 区 ,雨季 前 期 .雨季 和 雨季 后 期 
地 下 水 排泄 量 分 别 占 总 排泄 量 的 95. 1% .86. 1% 和 
75.2% 。 中 段 地 下 水 排泄 较 弱 。 雨 季 地 下 水 单 宽 排 
泄 量 明 显 降低 ,河水 中 基 流 比例 在 雨季 前 期 和 后 期 
分 别 为 68.9% 和 71.2% ,雨季 降 至 29.4% 。 
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Quantifications of spatial and temporal variations in groundwater 
discharge into a river using hydrochemical and isotopic tracers 
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Abstract: The Malian River, located in the southwestern portion of the Ordos Basin, Inner Mongolia, China, is a 
secondary tributary of the middle reaches of the Yellow River. It has experienced mounting pressure from exploration 
and utilization of water resources ,so the area has been facing serious water shortages and deterioration of water qual- 
ity. Finding accurate ways to quantify groundwater discharge to the river is critical for protecting and managing the 
water resources in the area. To study the spatial and temporal variations of groundwater recharge, we collected 75 
samples of surface water and groundwater in the stretch of the river downstream of the Malian River during pre-mon- 
soon , monsoon ,and post-monsoon seasons. A multi-tracer mass balance approach was applied to identify the loca- 
tion, amount, and seasonal variations of groundwater discharge. We used Cl concentrations and electrical conduc- 
tivity in conjunction with oxygen and hydrogen isotopes. We observed general decreases in Cl concentrations 
(387.20 -707.42 mg - L~') and EC (1 705.40 -3 123.60 ps - cm- ) ,as well as O ( -10.41%0 - —6. 30%) 
and 8D values ( —78.25%o - —54.88%c) as we moved downstream. The Cl” concentrations and EC values reached 
their minimum values during the monsoon season ,when 8'°O and 8D values peaked. All of the tracers had lower con- 
centrations or isotope ratios in the groundwater than they did in the river water. Groundwater values included Cl con- 
centrations (115.92 -255.96 mg - L~') ,EC (409.58 - 196.84 us + cm™') ,8"O ( - 11.19% — -9.11%o) , and 
òD ( -79.80%o - -68.28%o). No evident spatial or temporal variations in groundwater tracers were observed. All 
the groundwater samples had 8D and 8'*O values close to the Local Meteoric Water Line , indicating that groundwa- 
ter was recharged by precipitation and experienced limited evaporation. The river water was generally more enriched 
in 8D and O „which may have resulted from evaporation. Groundwater discharge rates varied spatially and tempo- 
rally. Groundwater in the upper and lower segments had much higher discharge rates ( accounting for 72.20% — 
95.07% of the total) than that in middle segment. The percentage of base flow in the river discharge was 68.89% 
during the pre-monsoon and 71.19% during the post-monsoon season. Groundwater inflows decreased sharply dur- 
ing monsoon season, with the ratio of base flow only accounting for 29.43% of the total river discharge. In general, 
the seasonal variations in groundwater discharge to the river were pronounced, while spatial variation was not evi- 
dent. The extent of interactions between river water and groundwater was mainly dominated by the geological and 
hydrogeological settings and conditions ,such as the aquifer thickness , water hydraulic gradient, and riverbed perme- 
ability ,so no spatial variations related to the stable conditions in the downstream direction were evident. Uncertain- 
ties remain in the model. However, when used correctly ,a multi-tracer method can be applied to provide detailed in- 
formation on the spatial distribution of groundwater discharge , especially in gaining rivers. This study could provide a 
theoretical basis and technical support for better understanding the interactions between river water and groundwa- 
ter,as well as the sustainable development of water resources and ecological and environmental protection. 
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